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FOREWORD 


The development of this computer program - referred to as TRAGEN - 
was supported under NASA Contract No. NASl-15497, by Langley Research 
Center, Hampton, Virginia. The project Technical Monitors were Samuel A. 
Morello, Kathy H. Samms, and Robert E. Shanks. At AMA, Inc. the project 
manager was John A. Sorensen, with engineering support provided by Mark H. 
Waters. Project programmers were Marianne N. Kidder, Quyen T.L. Nguyen, 
and Leda C. Patmore. 

This Users’ Guide describes the program input, program output, and 
general organization. Appendix A presents the technical material upon 
which the program is based. Appendix B presents a brief explanation of 
each of the program subroutines. 
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INTRODUCTION 


This doctunent is a technical description and a users' guide for a 
computer program - called TRAGEN - which is used to simulate an aircraft 
steered to follow a climbing, cruise, or descending profile or any 
sequential combination of these flight phases. Specifically, the pro- 
grcim simulates the longitudinal dynamics of a medium range twin- jet 
or tri-jet transport aircraft. For the climbing trajectory, the thrust 
is constrained to maximum value, and for descent, the thrust is set at 
idle. For cruise, the aircraft is held in the trim condition. 

For climb or descent, the aircraft is steered to follow either (a) 
a fixed profile which is input to the program or (b) a profile computed 
at the beginning of that segment of the run. For climb, the aircraft 
is steered to maintain the given airspeed as a function of altitude. 

For descent, the aircraft is steered to maintain the given altitude as a 
function of range-to-go. In both cases, the control variable is angle- 
of-attack. The given output trajectory is presented and compared with 
the input trajectory. Step climb is treated just as climb. 

For cruise, the Breguet equations are used to compute the fuel burned 
to achieve a given range and to connect given initial and final values of 
altitude and Mach number. 

TRAGEN is an acronym for trajectory gen eration. A companion program 
has been developed which generates optimum profiles which produce input 
to this program. This companion program - called OPTIM - uses optimiza- 
tion techniques and the energy state method to compute points on the pro- 
file. The conditions for climb and descent in OPTIM are consistent with 
those in TRAGEN. The users' guide for OPTIM appears as a separate document. 


Anon. , "OPTIM-Computer Program to Generate a Vertical Profile which 
Minimizes Aircraft Fuel Bum or Direct Operating Cost - Users' Guide", 
NASA CR-166061, March 1983. 
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The purposes of TRAGEN are three-fold: 


1. To verify the accuracy of near-optimum profiles generated by 
separate programs where simplifying assumptions are used to 
render the problem tractable. Specifically, TRAGEN can be 
used to verify the results of OPTIM. 

2. To compare the results of flying along a near-optimim profile 
with those generated by some other means. For example, the 
aircraft handbook specifies that the aircraft fly along 
profiles with fixed indicated airspeed and Mach numbers. 

3. To use as an evaluation tool for study of possible airborne 
implementation of autopilot/autothrottle flight management 
systems. For example, TRAGEN can be used to study the ability 
of the system to adapt to non-nominal flight conditions (e.g. , 
wind and atmospheric variations, change in destination). 

This users' guide is organized as follows: 

1) Section II presents, in concise form, the input cards and input 
data files that must be used to run the program. These are 
followed by a brief explanation of the options available to 
the user. 

2) Section III presents examples of the program output. This out- 
put consists of the profile followed by the aircraft and the 
reference profile. Without knowing the theory behind the con- 
struction of the program. Sections II and III enable the user 
to make runs and to interpret the results. 

3) Section IV presents the program layout in flowchart form. 

4) Appendix A gives a technical explanation of the aircraft 
equations of motion, how the steering control laws are 
generated, and the Breguet range equation. 

5) Appendix B describes the function of TRAGEN' s subroutines. 
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II 


INPUT DESCRIPTION 


TRAGEN is capable of running a multi-segment mission. Each segment 
requires up to five input cards and (optionally) two designated data files . 
The meanings of the variables on the input cards are given first. The 
program uses Unit 5 as the card input source. 

Card 1 

This card is the header that appears at the beginning of the segment. 
The input has an 8A10 format. 

Card 2 

This card consists of five integer variables used as flags to control 
the operation of the program, and one real environment variable. The 
input numbers are right- justified and have a 5I2,F10.0 format. They are: 

NSC JTRAJ lAC IWIND IPRT DTEMPK 

The meaning of each of these variables is as follows: 

NSC This is the mission segment control variable. The options 

presently available for NSC are: 

NSC = 2 : climb 

NSC = 3 : cruise 

NSC = 4 : descent 

NSC = 5 : end mission. 

This integer is read as the first item of input for each segment, and it 
is used to route the logic to the correct location within TRAGEN. For 
the initial segment, weight, airspeed, altitude, and time are all inputs, 
but they are all transferred within the program to match the end of a 
given segment with subsequent mission segments . 
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ITRAJ 


This flag is used to determine the source of the reference 
profile to be followed. Values are: 

ITRAJ =1 : Reference trajectory read in from a data file. 

ITRAJ =2 : Reference trajectory computed to follow a sequence 

of fixed Mach number and indicated airspeed segments. 

lAC This flag is used to select which aircraft model to use to 

generate the desired profile. Current values of lAC are: 

lAC = 2 : Medium-range three-engine jet transport aircraft, 

(nominal model) 

lAC = 3 : Medium-range two-engine jet transport aircraft. 

IWIND An arbitrary wind profile can be read in on Unit 7. It gives 
the simulated true wind speed and heading as a function of 
altitude. Options available for IWIND are; 

IWIND = 0 : No wind used. 

IWIND = 1 : Constant input wind profile used. 

IWIND = 2 : Separate wind profiles are used for climb, cruise 

and descent. 

NOTE: The profile is read only once for each mission (on the 

first leg) , although IWIND may be 0 or 1 for any mission 
segment. 

IPRT This flag is used to obtain additional printout of dynamic vari- 

ables during the integration process, as described in Section 
III, Output. Values are; 

IPRT =0 : No extra printout (normal mode) . 

IPRT = 1 : Reference trajectory printout included. 

IPRT = 2 : All printout included. 

DTEMPK This is the deviation from standard temperature, in degrees Celsius. 
Card 3 (Optional) 

This card has 3 real variables with format 3F10.3 and is read only 
for the initial segment. The variables are: 

HO VO WO 
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The meanings of these variables are as follows: 


HO 

VO 

WO 

Card 4 

This 

are: 

PSIG 

These are 
PSIG 

HP 

VF 

GRANGE 

VIAPl 

VIAP2 

RMP3 


Initial altitude, in ft. 

Initial Indicated airspeed, in kt. 
Initial aircraft weight, in lb. 


card has seven real variables with format 7F10.3. The variables 


HF VF GRANGE VIAPl VIAP2 PMP3 


defined as 


Aircraft heading over the ground in degrees. PSIG is used 
with the wind vector to compute the aircraft heading with 
respect to the air mass, in deg. 


Final altitude, in ft. 

Final Indicated airspeed, in kt. 

Cruise only: desired cruise range, mni. 


Ignored in cruise 
if reference 
trajectory is read in 


The desired Indicated airspeed in climbing to (or 
descending from) 10000 ft altitude, in kt. 

The desired Indicated airspeed in climbing from (or 
descending to) 10000 ft altitude to (from) intersection 
with Mach number RMP3. 

The desired Mach number in climbing from (or descending to) 
VIAP2 indicated airspeed to (from) Intersection with cruise 
altitude. 


required 
only for 
descent or 
initial 
climb and 
a reference 
trajectory 
is computed. 
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Card 5 - Climb and Descent only 

This card has seven real variables with format 7F10.3. The variables 

are: 

TSTOP DTI RKl RK2 RK3 RK4 ALFO 
The meanings of these variables are: 

TSTOP Time from the beginning of the integration to the stop 
time, in sec. 

DTI Integration step size, in sec. 

RKl Proportional gain used to convert airspeed error to 

angle-of-attack (6a) command, ("/ft/sec). 

RK2 Integral gain used to convert the integral of airspeed 

error to 6a command (°/ft). 

RK3 Proportional gain used to convert flight path angle error 

to 6a command (°/°). 

RK4 Integral gain used to convert the integral of flight path 

angle error to 6a command (°/° sec). 

ALFO The nominal value of angle-of-attack, in deg. 

Card 5 - Cruise Segment only. Optional 

This card has two real variables with format 3F10.0 and is read only 
if IWIND ^ 0 and the reference trajectory was not read in. This allows 
the use of changing winds for multiple cruise segments. The variables are: 

VWK PSIW 


Here, 



VWK 

Wind speed. 

in kt , 

PSIW 

Wind source 

direction. 
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Unit 7 - Wind Data (Optional) 

This data set is used when IWIND is set to 1 or 2. The input consists 
of the magnitude of the wind and the direction of its source as a function 
of altitude. The data format is (3F5.0 12). 

If IWIND = 1, a single wind profile applicable to the entire flight is 
read in. This profile consists of a set of n cards. Each card has four 
variables. 

HWIND(I) PSIW(I) VW(I) IE 

There is one card for each 1=1,2, — N, where N is the number of altitudes 
used for a given wind profile. The meanings of these variables are: 

HWIND(I) Beginning (lowest) altitude at which direction PSIW(I) and 
magnitude VW(I) apply. The program will interpolate for 
values of PSIW and VW when using altitudes between HWIND(I) 
and HWIND(I+1) . 

PSIW(I) Direction of the wind vector source in degrees (l.e., 270° 
represents a wind from the West) . 

VW(I) Magnitude of the wind vector, in kt. 

IE End-of-wind-table Indicator. If IE = 0 , the program will 

expect to read further wind data. If IE = 1 , the program 
assimes that a complete wind table has been read in. Note 
that when IE = 1, the corresponding altitude should be equal 
to or greater than any altitude the aircraft is expected to 
reach. 

If IWIND = 2, three wind profiles are read in, one each for climb, cruise, 
and descent (in that order). Each profile is as described under IWIND = 1. 
Each profile must end with a non-zero value for IE. If the reference 
trajectory is computed, the cruise wind is overwritten by data on Card 5. 
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Unit 11 Input reference data set 

These data are read in as the reference variables describing the 
nominal vertical profile followed by the aircraft during the climb and 
descent portions of the optimum trajectory. The data are obtained as 
output from Unit 11 in program OPTIM, although they could be obtained 
from any other source. 


The input consists of up to six binary records of the following form; 

Record 1; WORD, NWORD, 

WORD may be: CLIMB, CRUISE, or DESCEND. 

NWORD is the nimber of points stored for the 
specified flight segment. 

Record 2: An NWORD by 12 matrix of which only 10 columns 

are read. For example, for climb. Record 2 contains 
for JCLIMB = 1,..., NWORD, the following: 


CGRAF(J CLIMB, 1) = E 

CGRAF( JCLIMB, 2) = J 

CGRAF( JCLIMB, 3) = MACH 

CGRAF ( J CLIMB , 4 ) = VTASK 

CGRAF( JCLIMB ,5) = GAMMA 

CGRAF(JCLIMB,6) = FUELUZ 
CGRAF( JCLIMB, 7) = EPR 

CGRAF (JCLIMB, 8) = 0 

CGRAF( JCLIMB, 9) = TIME 

CGRAF(JCLIMB,10) = DIST 


Specific energy - ft 
altitude - ft 

Mach 

true airspeed - kt 

flight path angle - deg 

fuel burned - lb 

EPR settxng 

blank 

time 

range traveled - nmi 


The same variables are stored in SGRAF for the cruise and DGRAF for the 
descent portion. The descent profile is generated backwards in time in 
OPTIM. 
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Ill 


OUTPUT DESCRIPTION 

The output of TRAGEN is compact and generally self explanatory. The 
output is printed using Unit 6. The input quantity IPRT controls the amount 
of output. 

The first output for each mission segment consists of printing the 
input. A typical example for a climb segment is shown in Table la. An 
example for a cruise segment is shown in Table lb. The first line is the 
header which is used to identify the run. The next several lines print 
out and explain the run control flags. The last lines print out the 
real number program control variables. Definitions of these flags and 
variables are the same as presented in Chapter II. 

The cruise segment example also contains the initial and final per- 
formance values. The initial values, in this case, have been constructed 
from the results of the previous climb or cruise segment. 

The next set of output is dependent upon whether the reference pro- 
files is read in or computed. If ITRAJ =2, the profile is computed. 

If IPRT = 2, printout is Included which indicates consecutive variables 
of this computation. An example of this output is shown as Table 2. 

The variables that are printed out are: 

VIAS indicated airspeed (ft/sec) 

RM Mach number 

D drag (lb) 

TH thrust (lb) 

W weight (lb) 

EDT energy rate (ft/sec) 

FDT fuel flow rate (Ib/hr) 

T time (sec) 

R range (ft) 

H altitude (ft) 

V true airspeed (ft/sec) 
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Table la. Input Data Printout, Climb 



INSC : 3 

CRUISE 

PHASE 

0 


CRUISE SEGMENT *1 

IWIND 

IPRT 

lAC 

1 

0 

3 

OIUlIND : 1 

INPUT 

UIND 


OIPRT : 0 NO EXTRA PRINT 

OTAC : 3 TUO-JET AIRCRAFT HODEL 

DESIRED AIRCRAl-r HEADING (DEG FROM NORTH) 
WIND SPEED (KD 

WIND DIRECTION (SOURCE HEADING) 
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CRUISE SEGMENT PERFORMANCE// INPUT DATA 


INITIAL CONDITIONS 


Al TI rUDF 

(FT) 

31603. 

MACH NO 


.760 

WEIGH I 

(LBS) 

92489. 

TIME 

(HRS) 

. 231 

RAN(3E 

(NMI ) 

88. 

FINAL CONDEIIOHS 

ALinUDE 

(FI ) 

32000. 

MACH NO 


.760 

SEGMENT RANGE 

(NMI ) 

250. 


oumn DATA 


AlIITUDE (FD 

MACH NO 

UEIUIIT a.BS; 

rrME (HRS) 

RANGE tNMI) 

RANGE FACrOR (NMH 


m-OFN FND 

31603. 3/000. 

./(f.0 . /60 

y'14 80, GV23V. 

.231 .B6S 

88. 338. 

/899. 7800. 


Table lb. Input Data Printout, Cruise 
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VlAS* Rri# 0* TH« U« E0T« EOT 421.993 .441 9310.742 23636.067 146877.597 49.760 17210.929 

T> R. H, V* GAH» E» ft EP* VW 216.70 99937.87 9900.00 475.97 9.11 13011.90 1130.01 1.89 0.00 

VIAS# RH> 0* TH» M* EOT« POT 421.993 .445 9306.590 23327.584 148819.992 45.138 16993.425 

It H* V« GAN* E* F* EP« VM 226.77 101298.59 10000.00 479.10 5.00 13564.16 1237.70 1.90 0.00 

VIAS* RN* 0* TH# U« EOT* FOT 906.343 .928 10437.790 22974.488 148762.296 47.966 17542.194 

T* R* H* V* GAN* E» ft EP* VW 261.25 118323.08 10000.00 569.17 0.00 15030.30 1391.03 1.90 0.00 

VIAS* RN* 0* TH* W* EOT* FOT 906.343 .533 10429.600 22694.833 148608.971 47.331 17330.996 

T* R* H* V* GAN* E* F* EP* VW 273.20 125150.46 10900.00 573.30 4.22 19603.56 1449.27 1.87 0.00 

VIAS* RN* 0* TH* W* EOT* FOT 506.343 .938 10416.914 22417.644 143550.731 46.651 17122.726 

T* R* H* V* GAN* E* F* EP* VW 285.34 132134.61 11000.00 977.47 4.12 16176.10 1907.70 1.88 0.00 

VIAS* RN* 0* TH* W* EOT* FOT 506.343 .543 10407.893 22141.652 148492.295 49.964 16916.367 

T* R* H* V* GAN* E* F* EP* VW 297.68 139287.74 11500.00 581.66 4.02 16753.87 1566.41 1.88 0.00 

VIAS* RN* 0* TH* W* EOT* FOT 906.343 .548 10398.732 21866.092 148433.593 45.269 16711.933 

T* R* H* V* GAN* E* ft EP* VW 310.24 146617.00 12000.00 985.93 3.93 17330.91 1625.40 1.89 0.00 

VIAS* RN* 0* TH* W* EOT* FOT 906.343 .953 10389.430 21593.398 146374.600 44.568 16509.433 

T* R* H* V* GAN* E* F* EP* VW 323.01 154129.97 12500.00 590.22 3.83 17909.26 1684.71 1.89 0.00 

VIAS* RN* 0* TH* W* EOT* FOT 506.343 .558 10379.962 21320.850 148315.293 43.858 16308.875 

T* R* H* V* GAN* E> F* EP* VW 336.02 161834.65 13000.00 594.59 3.73 18488.92 1744.35 1.90 0.00 

VIAS* RN* 0* TH* W* EOT* FOT 906.343 .563 10370.386 21049.372 148255.647 43.141 16110.267 

T* R* H* V* GAN* E* ft EP* VW 349.27 169739.81 13500.00 598.92 3.64 19069.93 1804.37 1.90 0.00 

VIAS* RN* 0* TH* W* EOT* FOT 506.343 .569 10360.639 20779.186 148195.633 42.416 15913.616 

T* R* H* V* GAN* E* ft EP* VW 362.77 177854.71 14000.00 603.33 3.99 19652.31 1864.78 1.91 0.00 

VIAS* RN* 0* TH* W* EOT* FOT 906.343 .574 10350.737 20510.327 148135.222 41.684 15718.926 

T* R* H* V* GAN* E* F* EP* VW 376.53 186189.04 14500.00 607.79 3.49 20236.08 1925.62 1.91 0.00 

VIAS* RN* 0* TH* W* EOT* FOT 506.343 .579 10340.677 20242.826 148074.383 40.945 15526.204 

K. T* R* H* V* GAN* E* ft EP* VW 390.57 194753.12 15000.00 612.26 3.36 20821.27 1986.92 1.92 0.00 

N3 VIAS* RN* D» TH* W* EOT* FOT 506.343 .989 10330.496 19976.717 148013.085 40.199 15339.452 

T* R* H* V* GAN* E* ft EP* VW 404.90 203557.95 15500.00 616.82 3.27 21407.90 2048.71 1.92 0.00 

VIAS* RN* 0* TH* W* EOT* FDT 506.343 .590 10320.070 19712.029 147951.294 39.447 15146.675 

T* R* H* V* GAN* E* ft EP* VW 419.53 212619.28 16000.00 621.40 3.18 21995.99 2111.03 1.93 0.00 

VIAS* RN* 0* TH* W* EOT* FOT 906.343 .996 10309.915 19448.793 147688.974 38.687 14959.874 

T* R* H* V* GAN* E* F* EP* VW 434.47 221937.64 16500.00 626.03 3.09 22585.58 2173.91 1.93 0.00 

VIAS* RN* 0* TH* W* EOT* FOT 906.343 .601 10298.769 19189.328 147826.088 37.931 14773.852 

T* R* H» V* GAN* E* F* EP* VW 449.75 231530.46 17000.00 630.70 3.00 23176.69 2237.40 1.94 0.00 

VIAS* RN* 0* TH* W* EOT* FOT 906.343 .607 10287.887 18942.418 147762.595 37.216 14587.731 

T* R* Hf V* GAN* E* F* EP* VW 465.38 241429.56 17500.00 635.41 2.91 23769.35 2301.52 1.94 0.00 

VIAS* RN* 0* TH* W* EOT* FOT 506.343 .613 10276.807 16696.704 147698.475 36.494 14403.698 

T* R* H* V* GAN* E* ft EP* VW 481.34 251613.04 18000.00 640.17 2.82 24363.57 2366.22 1.95 0.00 

VIAS* RN* 0* TH* W* EOT* FOT 506.343 .618 10265.547 18452.216 147633.775 35.765 14221.750 

T* R* H* V* GAN* E* ft EP* VW 497.67 262104.01 18500.00 644.97 2.74 24959.40 2431.99 1.95 0.00 

VIAS* RN* 0* TH* W* EOT* FOT 906.343 .624 10254.105 18208.963 147566.452 39.029 14041.884 

T* R* H* V* GAN* E* ft EP* VW 514.38 272918.67 19000.00 649.82 2.66 29556.86 2497.54 1.96 0.00 

VIAS* RN* 0* TH* Wf EOT* FOT 906.343 .630 10242.479 17967.034 147502.459 34.286 13864.095 


Table 2. Computed Reference Trajectory Variables 



GAM flight path angle (deg) 

E specific energy (ft) 

F fuel burned (lb) 

EP EPR setting 

VW wind speed (ft/sec) 

These data are useful for checking over the computation of the reference 
profile. 

Next, if IPRT = 1 or 2, the reference profile is printed. This is as 
computed by the TRAGEN subroutine REFCOM or as read in from Unit 11. An 
example of an input reference trajectory Is shown as Table 3. 

Next, the vertical wind profile is printed as shown in Table 4. This 
is used when there is a non-zero wind, and the flag IWIND is set to 1. 

Next, the wind shear data is printed as shown in Table 5. The wind 
data are taken from Table 4 and used to compute the North and East com- 
ponents of shear every 2000 ft. These shear data are used by the pro- 
gram to compute the real-world (or actual) aircraft longitudinal component 
of wind as a function of altitude. This simulated actual shear may differ 
from what is recorded on the reference trajectory data set. 

Finally, the main results of the program are printed as shown in 
Tables 6 and 7. Table 6 shows the first page of printout for a descending 
profile generated by using the autopilot logic of STEER2 (described later) . 
This print consists of a sequence of two lines of variables shown at nearly 
equal time points. The first line of variables are those generated by 
integrating the aircraft equations of motion (described in Appendix A) . 

The second line contains the same variables as obtained from the input 
nominal reference trajectory. The variables are the same as those printed 
for the nominal reference trajectory. 

This print enables the user to directly compare the performance obtained 
from the integrated equations of motion with that predicted from another 
source. (This may be the optimum performance as computed by program OPTIM) . 
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descend OPTIHIZATION reference TRAJECTORr DATA 


TINE 

RANGE 

altitude 

AIRSPEED 

GANNA 

energy 

FUEL 

EPR 

wind V 

SEC 

FT 

FT 

FT/S 

OEG 

FT 

LB 


FT/S 

•87A.6A6 

-666602.110 

33118.839 

736.365 

-3.256 

61538.163 

0.000 

1.803 

-56.836 

-869.555 

-661161.260 

32922.101 

736.266 

-3.890 

61339.163 

1.782 

1.303 

-55.911 

-86A.487 

-657706.697 

32688.562 

735.626 

-6,075 

61089.163 

3.555 

1.803 

-57,171 

-859. A39 

-656296.311 

32665.662 

735.225 

-6,105 

60839.163 

5.323 

1.803 

-58,662 

-85A.A16 

-650907.939 

32202.^02 

736.920 

-7,661 

60589,163 

7.080 

1.803 

-59,727 

-869. AZO 

-667567.153 

31763.686 

763.151 

-3.960 

60339.163 

8.829 

1.803 

-61.556 

-866.599 

-666273.069 

31538.009 

762.088 

-6,129 

60089.163 

10.517 

1.803 

-62.300 

-839.791 

-661016.688 

31302.962 

761.1639 

-6.135 

39839.163 

12.199 

1.803 

-63.085 

-835.095 

-637781,890 

31069.108 

760..737 

-6.205 

39589.163 

13.876 

1.803 

-63.851 

-830.238 

-636567,697 

30832.790 

760.062 

-5.538 

39339.163 

15.563 

1.803 

-66.617 

-825.696 

-631373.696 

30523.103 

762.736 

-6,581 

39089.163 

17.203 

1.803 

-65.610 

-820.827 

-628225.585 

30199.913 

767.625 

-6,726 

38839.163 

18.837 

1.803 

-66,757 

-816.276 

-629138.800 

29795.970 

752.517 

-8,716 

38589.163 

20,630 

1.803 

-67.627 

-811.838 

-622138.356 

29336.073 

/ 750.799 

-6.668 

38089.163 

21.983 

1.803 

-67.816 

-803.026 

-616165.977 

28867.867 

769.1635 

-6,517 

37589.163 

25.067 

1.803 

-68.207 

-796.290 

-610220.636 

28399.776 

763.062 

-3.908 

37089,163 

28.125 

1.803 

-68.601 

-785.629 

-606369.597 

28000.061 

763.733 

-3.755 

36589.163 

31.156 

1.303 

-68,937 

-776.965 

-398566.009 

27617.815 

738.617 

-6.339 

36089,163 

36.195 

1.803 

-69,258 

-768.212 

-392729.897 

27176.669 

736.055 

-5.780 

35589.163 

37.273 

1.803 

-69.629 

-759.515 

-386969,673 

26991.356 

739.769 

-6,178 

35089.163 

60.358 

1.803 

-70.121 

-751.063 

-381306.956 

25980.321 

766.587 

-6.360 

36589.163 

63.387 

1.803 

-70.668 

-762.889 

-375816.667 

29368.566 

769,1605 

-6,555 

36009.163 

66,353 

1.803 

-72.211 

-736.986 

-370685.589 

26756.029 

756.226 

-6,768 

33589.163 

69.261 

1.803 

-73.756 

-727.333 

-365302.812 

26162.770 

759.066 

-6,915 

33089.163 

52.162 

1.803 

-75,302 

-719.915 

-360260.038 

23528.767 

763.365 

-6.929 

32589.163 

55.039 

1.803 

-72.693 

-712.727 

-359289,522 

23101.833 

760,779 

-6.979 

32089.163 

57.967 

1.803 

-59.622 

-705.563 

-350363.117 

22670.937 

757.368 

-6.981 

31569.163 

60.906 

1.803 

-56.726 

-698.367 

-369602.738 

22260.366 

756.893 

-2,817 

31089.163 

63.910 

1.803 

-63.828 

-691.199 

-360688.766 

21998.580. 

763.t797 

-3.056 

30589.163 

66.965 

1.803 

-62.206 

-683.801 

-335687.522 

21731.725 

733.636 

-5.181 

30089.163 

70.119 

1.803 

-61.262 

-676.192 

-330392.263 

21269.698 

731.966 

-5,162 

29589,163 

73.368 

1.803 

-59,559 

-668.627 

-325326.787 

20812.099 

730.098 

-5,236 

29089.163 

76.655 

1.8)3 

-57.875 

-661.100 

-320287.616 

20350.670 

728.1606 

-6,593 

20589,163 

79.980 

1.803 

-56.161 

-653.615 

-315282.869 

19968.393 

726.062 

-5.682 

28089.163 

63.363 

1.803 

-56.722 

-666.092 

-310267.783 

19669.600 

726.017 

-5.762 

27589,163 

86.763 

1.803 

-53.617 

-638.663 

-309305.758 

18950.613 

723.972 

-5.706 

27089,163 

90.219 

1.803 

-52.028 

-631.326 

-300396.666 

18659,881 

723.551 

-5.368 

26589.163 

93.710 

1.803 

-50.585 


Table 3. The Input or Computed Reference Trajectory 


WIND DATA 


FTI, 

VW(KN0TS)> 

VW(FT/SECI# 

PSIW(OEG) 

0. 

0*00 

0*00 

205* 

2000. 

4*00 

6*75 

205* 

A03D. 

9.00 

15*19 

205. 

600C« 

14*00 

23*63 

205* 

803C. 

13*00 

21*94 

240, 

lOOOC. 

14.00 

23*63 

275* 

12000* 

23*00 

38*62 

285* 

14030* 

22*00 

37*13 

275, 

16000. 

27*00 

45*57 

270* 

1800C. 

31.00 

52*32 

290* 

2000C* 

33*00 

55.70 

289* 

2200C. 

37*00 

62*45 

275* 

2400C* 

45*00 

75*95 

275* 

26000* 

42.00 

70*89 

275* 

2800C* 

41.00 

69*20 

275* 

30000* 

40*00 

67,51 

275* 

32000* 

36*00 

60*76 

270* 

3400r* 

30*00 

50*63 

280. 

36000* 

28*00 

47*26 

300* 

3800C* 

31*00 

52*32 

310. 

40000* 

31*00 

52*32 

310* 

42000* 

31*00 

52*32 

310* 

44000* 

31*00 

52*32 

310,; 

4600C* 

31*00 

52*32 

310* 


Table A . Vertical Wind Profile 
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WIND SHEAR DATA 


ALT 

VU 

PSjW 

D(WX)/DH 

0{wy) /OH 

C. 

0.00 

205. 

-.003059 

-.001427 

2000. 

6.75 

205. 

-.003824 

-.001733 

4000, 

15.19 

205 • 

-.003824 

-.001733 

600C . 

23.63 

205. 

.005222 

-.00.4503 

800C. 

21.94 

240. 

.006515 

-.002269 

lOOOC. 

23.63 

275. 

.003994 

-.006979 

1200C. 

33.62 

285. 

-.003406 

.000253 

1400C. 

37.13 

275, 

-.001618 

-.004290 

1600D. 

45.57 

270. 

.008948 

-.001798 

leooc. 

52.32 

290. 

-.004112 

-.002842 

20000. 

55.70 

280. 

-.002115 

-.003630 

2200C. 

62.45 

275. 

.000583 

-.006726 

2400C. 

75.95 

275. 

-.000221 

.002522 

2600C. 

70.89 

275. 

-.000074 

•000841 

2800C. 

69.20 

275, 

-.000074 

.000341 

3000C. 

67.51 

275. 

-.002942 

•003247 

3200C. 

60.76 

270, 

•004396 

.005443 

3400C. 

50.63 

28C. 

.007413 

.004469 

3600C. 

47.26 

300. 

.005001 

•900423 

3800D. 

52.32 

310. 

0.000000 

9.000000 

4000C. 

52.32 

310. 

0.000000 

0.000000. 

4200C. 

52,32 

310. 

0.000000 

0.000000 

4400C. 

52.32 

310. 

0.000000. 

0.000000 


Table 5, Wind Shear Data as Functions of Altitude 
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1 CHM6 TRAJECTORY COMPARISON USINC CONTROL OPTION 1 


0 

TIME 

RANGE 

ALTITUDE 

AIRSFECP 

GAMMA 

ENERGY 

FUEL 

EPR 

HIND 

MACH 

ALPHA 


SEC 

NMI 

FT 

FT/S 

KT 

DEG 

FT 

LB 


FT/S 

NO 

DEG 

OACT 

0.000 

0.000 

5.0 

354.440 

210.000 

0 000 

1950.7 

0.000 

1 . 849 

0.000 

.317 

4.600 

REF 

0.000 

0.000 

0.0 

354.440 

210.000 

0.000 

1950.7 

0.000 

1.849 

0.000 

.317 


OACT 

8.500 

.527 

-8.9 

397.992 

235.804 

.975 

2452.7 

35.817 

1.840 

0.000 

.356 

4.669 

REF 

8.203 

410 

2.9 

397.042 

235.241 

054 

2450.7 

34.730 

1.840 

0.000 

.356 


OACT 

18.500 

1 . 209 

192.1 

430.243 

254.912 

3 335 

3068.7 

78.269 

1.835 

0.000 

.386 

3.962 

REF 

16.145 

.744 

188.6 

421.763 

249.888 

3.274 

2950.7 

68.409 

1.837 

0.000 

.378 


OACT 

28.500 

1.924 

677.9 

440. 101 

260.752 

8. 355 

3688.0 

120.631 

1.838 

0.000 

.395 

3.857 

REF 

24.181 

1.499 

672.0 

423.026 

250.636 

8.188 

3450.7 

102.079 

1.842 

0.000 

.380 


OACT 

35 500 

2.427 

1091.8 

442.751 

262.323 

6.749 

4138.2 

149.937 

1.842 

0.000 

.398 

3.394 

REF 

32.255 

2.059 

1098.5 

428.586 

253.930 

7 126 

3950.7 

135. 585 

1.846 

0.000 

.385 


OACT 

44.000 

3.044 

1572.7 

446.398 

264.484 

7.757 

4669.5 

185. 171 

1.847 

0.000 

.402 

3 485 

REF 

40.409 

2.631 

1563.7 

431.194 

255.475 

7.626 

4450.7 

169.042 

1.851 

0.000 

.388 


OACT 

51.500 

3.591 

2011.3 

448.596 

265.786 

7.344 

5138 6 

215.901 

1.852 

0.000 

.405 

,3 . 523 

REF 

48.A38 

3.211 

2025.4 

434.039 

257. 161 

7.453 

4950.7 

202.440 

1.856 

0.000 

.391 


OACT 

59.500 

4. 178 

2471.0 

451.051 

267.241 

7.294 

5632. 7 

248.329 

1.858 

0.000 

.407 

3.516 

REF 

56.935 

3.801 

2486.9 

436.891 

258.851 

7.337 

5450.7 

235.775 

1.862 

0.000 

.395 


OACT 

68.000 

4.805 

2954.3 

453.681 

268.799 

7.163 

6153.0 

282. 421 

1.864 

0.000 

.411 

3.538 

REF 

65.303 

4.399 

2947.7 

439.766 

260.554 

7.218 

5950 7 

269.052 

1 . 868 

0.000 

.398 


OACT 

76.000 

5. 399 

3404.6 

456.157 

270.265 

7.049 

6618.2 

314.176 

1.870 

0.000 

.413 

3.550 

REF 

73.743 

5.007 

3408.0 

442.663 

262.271 

7.100 

6450.7 

302.274 

1.874 

0.000 

.401 


OACT 

84.500 

6.034 

3877.8 

458.801 

271.832 

6.931 

7149.0 

347.568 

1.876 

0.000 

.417 

3.563 

REF 

82.257 

5 . 624 

3867.8 

445.503 

264.001 

6.983 

6950.7 

335. 445 

1.380 

0.000 

.405 



Table 6 


Actual and Reference Optimum Descent Profiles 



Table 7 shows the first page of a more dense printout for a climb 
profile, obtained by setting the flag IPRT =2. The print contains the 
same two lines of variables as shown in Table 6 at various reference 
points in time. It also shows three different types of secondary print- 
out between these reference 'points. From Table 7, these are: 

Type 1: DGDH DVDH GCl VCl DT ICFL 

These are variables generated by the command logic (STEERl) to 
produce steering commands over the next period of time. For 
the climb profile governed by STEERl, these variables are; 

DGDH the gradient of flight path angle with altitude (deg/ft), 

DVDH the gradient of airspeed with altitude (ft/sec/ft), 

GCl constant term in the flight path-angle command (deg) 

VCl constant term in the airspeed command (ft/sec), 

DT time Interval over which the command applies (sec) , 

ICFL flag used to command (a) airspeed and flight path angle 

(ICFL = 1) , or (b) only flight path angle (ICFL = 2) . 

Type 2: VCM GCM DVD DGD DALF ALF VEI GEI 

These variables are generated four times per integration step. 
Their meanings are; 

VCM commanded airspeed (ft/sec), 

GCM commanded flight path angle (deg) , 

DVD error in airspeed (deg) , 

DGD error in flight path angle (deg) , 

DALF commanded incremental angle-of-attack 6a (deg) , 

ALF total commanded angle-of-attack (deg) , 

VEI Integral of airspeed error (ft) , 

GEI integral of flight-path-angle error (deg-sec) . 
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CLIMB TRAJECTORY COMPARISON USING CONTROL OPTION 1 

TIME RANGE ALTITUDE AIRSPEED GAMMA ENERGY FUEL EPR UIND M ALPHA 

SEC FT FT FT/S DEG FT LB FT/S DEG 

act 0.000 0.000 5.000 350.154 0.000 1903.843 0.000 1.850 0.000 

ref 0.000 0.000 0.000 350.154 0.000 1903.843 0.000 1.850 0.000 

DGDHfDUDHfGClfVClfDTiICFL 0.0000 10.495A 0.0000 387.5684 100.0000 2 

UCMf GCMt DVDf DGD» DALF. ALFf UEI* GEI 0.00 0.00 0.00 0.00 0.00 4.00 0.00 0.00 

L.UfTH.D.MAS.GAM.ALF 69266.1 95000.0 21555.5 4201.2 2952.7 0.0 .1 

HDD*UAD.XD»HDT.HD.VAiXfH -8.2 5.9 350.2 0.0 0.0 350.2 0.0 5.0 

VCM» GCMf DUDj DGDf DALFf ALF» VEIf GEI 0.00 0.00 0.00 .20 .08 4.08 0.00 0.00 

LfM»TH»DfMAS»GAMiALF 70763.9 94999.4 21548.9 4269.5 2952.7 -.0 .1 

HDDfUADfXD.HDT.HDrWA.XfH -7.7 5.9 351.0 -1.2 -1.2 351.0 52.5 5.0 

VCMf GCMj DUDr DGDf DALFf ALFf UEIf GEI 0.00 0.00 0.00 .44 .21 4.21 0.00 .14 

LfUfTHfDfMASfGAMf ALF 73122.0 94998.4 21539.1 4378.6 2952.6 -.0 .1 

HDDfVADfXDfHDTfHDfWAfXfH -6.9 6.0 352.4 -2.7 -2.7 352.4 131.9 4.1 

UCMf GCMf DVDf DGDf DALFf ALFf VEI f GEI 0.00 0.00 0.00 .66 .26 4.26 0.00 -.04 

LfUtTHfDfMASfGAMfALF 74126.6 94997.9 21533.5 4426.3 2952.6 -.0 .1 

HDDfUADfXDfHDTfHDfUAfXfH -6.6 6.1 353.0 -4.1 -4.1 353.1 175.1 5.2 

Mf Af EPf THf GAf LfD .3162 4.2562 1.8494 21533.5410 -.6609 74126.6269 4426.2809 

AZf AXf VXf VZf UDf Vf Xf H -6.5957 6.0673 353.0401 -4.0725 -4.2034 353.1304 175.8139 

UCMf GCMf DVDf DGDf DALFf ALFf VEIf GEI 0.00 0.00 0.00 .60 .27 4.27 0.00 .16 

LfUfTHfDfMASfQAMfALF 74358.5 94997.9 21533.7 4436.8 2952.6 -.0 .1 

HDDfVADfXDfHDTfHDfVAfXfH -6.5 6.0 353.1 -3.7 -3.7 353.1 175.8 4.0 

VCMf GCMf DVDf DGDf DALFf ALFf VEIf GEI 0.00 0.00 0.00 .75 .35 4.35 0.00 .25 

LfUfTHfDfMASfGAMf ALF 75911.3 94997.3 21527.2 4510.5 2952.6 -.0 .1 

HDDfVADfXDfHDTfHDfVAfXfH -6.0 6.1 354.0 -4.7 -4.7 354.0 228.8 3.5 

VCMf GCMf DVDf DGDf DALFf ALFf VEIf GEI 0.00 0.00 0.00 .93 .47 4.47 0.00 .49 

LfUfTHfDfMASfGAMfALF 78245.6 94996.3 21517.6 4623.0 2952.6 -.0 .1 

HDDfVADfXDfHDTfHDfVAfXfH -5.2 6.1 355.4 -5.8 -5.8 355.4 308.9 2.0 

^ VCMf GCMf DVDf DGDf DALFf ALFf VEIf DEI 0.00 0.00 0.00 1.12 .53 4.53 0.00 .41 

LfUfTHfDfMASfGAMfALF 79427.5 94995.8 21512.2 4680.7 2952.6 -.0 .1 

HDDfVADfXDfHDTfHDfVAfXfH -4.8 6.2 356.1 -6.9 -6.9 356.1 352.4 2.5 

Mf Af EPf THf GAf LfD .3190 4.5285 1.8488 21512.1759 -1.1151 79427.5039 4680.7179 

AZf AXf VXf VZf UDf Vf Xf H -4.8132 6.1717 356.0615 -6.9309 -4.2069 356.1770 353.1193 1.462 

VCMf GCMf DVDf DGDf DALFf ALFf VEIf GEI 0.00 0.00 0.00 1.05 .53 4.53 0.00 .57 

LfUfTHfDfMASfGAMfALF 79516.2 94995.8 21512.3 4685.1 2952.6 -.0 .1 

HDDfVADfXDfHDTfHDfVAfXfH -4.8 6.1 356.1 -6.5 -6.5 356.2 353.1 1.5 

VCMf GCMf DVDf DGDf DALFf ALFf VEIf GEI 0.00 0.00 0.00 1.16 .61 4.61 0.00 .73 

LfUfTHfDfMASfGAMfALF 81055.6 94995.2 21506.0 4761.0 2952.5 -.0 .1 

HDDfVADfXDfHDTfHDfVAfXfH -4.3 6.2 357.0 -7.2 -7.2 357.1 406.5 .5 

VCMf GCMf DVDf DGDf DALFf ALFf VEIf GEI 0.00 0.00 0.00 1.27 .72 4.72 0.00 1.04 

LfUfTHfDfMASfGAMfALF 83271.4 94994.2 21496.0 4872.2 2952.5 -.0 .1 

HDDfVADfXDfHDTfHDfVAfXfH -3.5 6.2 358.4 -7.9 -7.9 358.5 487.3 -1.5 

VCMf GCMf DVDf DGDf DALFf ALFf VEIf GEI 0.00 0.00 0.00 1.42 .78 4.78 0.00 1.05 

LfUfTHfDfMASfGAMfALF 84563.6 94993.7 21490.6 4937.3 2952.5 -.0 .1 

HDDfVADfXDfHDTfHDfVAfXfH -3.1 6.2 359.1 -8.9 -8.9 359.2 531.2 -1.5 

Mf Af EPf THf GAf LfD .3218 4.7793 1.8481 21490.6126 -1.4210 84563.6122 4937.2793 

AZf AXf VXf VZf UDf Vf Xf H -3.0736 6.2089 359.1187 -8.9086 -4.2105 359.2590 531.9374 -2.314 

VCMf GCMf DVDf DGDf DALFf ALFf VEIf GEI 0.00 0.00 0.00 1.35 .78 4.78 0.00 1.18 


4.000 


4.044 


Table 7 


Climb Profile With Secondary Printout 



Type 3: M A EP TH GA L D 

AZ AX VX VZ WD V X H 

These variables indicate the state of the aircraft at the end 
of each integration step. The meanings of these variables are: 

M Mach number, 

A angle-of-attack (deg), 

EP EPR setting , 

TH thrust (lb) , 

GA flight path angle (deg) , 

L lift (Ib)^ 

D drag (lb), 

2 

AZ vertical acceleration (ft/sec ), 

2 

AX horizontal acceleration (ft/sec )^ 

VX ground speed (ft/sec), 

VZ altitude rate (ft/sec), 

WD fuel burn rate (Ib/sec), 

V airspeed (ft/sec), 

X distance-to-go or range (ft), 

H altitude (ft). 

If IPRT = 2, secondary output is also created for the descent trajectory. 
Again, there are three types of secondary printout. These are: 

Type 1: HP H RP X DHDX GCl DT 

These are variables used and generated by the command logic 
(STEER2) to produce steering commands over the next period 
of time. For the descent profile governed by STEER2, these 
variables are: 

HP the next reference altitude point (ft)^ 

H current measured altitude (ft) , 
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RP 

the next reference range point (ft) , 



X 

current measured range (ft). 



DHDX 

computed gradient of dh/dx. 



GCl 

commanded inertial flight path angle 

(deg) , 


DT 

time interval over which the command 
case, the command GCL is regenerated 
than RP. 

applies (sec) , 
when X becomes 

In this 
greater 


Type 2: GCM GMG DALF ALF GEI 

These variables are generated four times per integration step. 
They are the result of a fourth-order Runge-Kutta-Gill integra- 
tion method which is subroutine GO. The meanings of these 
variables are: 

GCM commanded Inertial flight path angle (deg) , 

GMG actual inertial flight path angle (deg) , 

DALF commanded Incremental angle-of-attack 6a (deg) , 

ALF total commanded angle-of-attack (a^ + 6a) (deg) , 

GEI Integral of flight-path-angle error (deg-sec) . 

Type 3 is the same format as produced for climb. 
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IV 


PROGRAM ORGANIZATION AND SUBROUTINES 

This section gives a brief overview of the process used in TRAGEN 
to generate steering conmands and to integrate the equations of motion 
for an aircraft following a given reference trajectory. This section 
also contains a brief description of each of TRAGEN 's thirty-nine sub- 
routines. The technical details upon which the program is based are 
presented in Appendix A. A more detailed description of the subroutines 
is presented in Appendix B. 

Figure 1 is a flow chart of the steps followed by TRAGEN to simulate 
an aircraft steered to follow an input or computed reference profile. The 
steps followed by the program are as follows: 

1. Read in the control flags, reference trajectory, prevailing 
wind model, and program control parameters. Use these data 
to initialize the program variables. 

2. If the desired segment is a cruise, calculate the cruise 
performance. Then return to step 1. 

3. If mission segment is a climb, and not the initial segment, 
perform a step climb subsegment. 

4. If a reference trajectory is to be computed (ITRAJ equals 2), 
compute a reference climb or descent trajectory. This consists 
of incrementing the altitude in steps of 500 ft, and computing 
the associated aircraft variables so that the desired true air- 
speed is maintained. The desired true airspeed is computed 
from the input sequence of Indicated airspeeds VIAPl and VIAP2 
and the Mach number RM3. This speed profile is similar to those 
specified in a typical pilot's handbook. 

5. Start the simulation update process. This consists of first 
writing the state variables as determined from integrating the 
aircraft equations of motion. These are written along with 
similar variables taken from the input reference profile. The 
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1 Read in 

• Control Flags 
Reference Trajectory 
Wind Data 

• Aircraft Data 



Figure 1. Macro Flow Chart Trajectory Generation Program 
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reference profile consists of a sequence of discrete points 
where the specific energy or the altitude changes in steps of 
500 ft. The steering algorithm commands a continuous trajectory 
which connects these points. After the actual and reference 
trajectory states are written, the program determines whether 
the stopping time (TSTOP) has been reached. 

6. If the simulation is to continue, the program next calls the 
appropriate steering algorithm. One method is available to gener- 
ate angle-of-attack commands to maintain airspeed as a function 
of altitude for the climb profile. A different method is avail- 
able to generate angle-of-attack commands to maintain altitude 
and flight path angle as a function of range-to-go for the descent 
profile. This block is called when the aircraft is near each 
reference point along the input profile. 

7. This step Integrates the equations of motion of the aircraft. 

The equation derivatives are computed four times per Runge 
Kutta integration step (Input DTI). Integration continues until 
the simulated aircraft reaches the next reference value of 
altitude or airspeed during climb or the next reference value 

of range during descent. At the end of the segment (next 
reference point reached), the program loops back to Step 3. 

8. Return to step 1. 

These steps are explained in more detail in Appendix A. 

TRAGEN is programmed in FORTRAN, and it consists of the main executive 
program and thirty-nine subroutines and functions. These thirty-nine sub- 
routines are called to execute the steps depicted in Fig. 1. Explanations 
of the program and its subroutines are presented in Appendix B. 

The program subroutines can be grouped into four categories: 

1. models of airborne software used to compute and generate 
steering commands, 

2. aerodynamic, propulsion, and flight dynamics models, 

3. flight condition models, and 

4. utility programs. 
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Under Category 1, the subroutines are: 


ACRUSE 

CMACTL 

CREWZ 

REFCOM 

SETREF 

STEERl 

STEER2 

VTCM 


Sets up program variables to initialize cruise segment. 

Initializes variables and controls computational flow 
for Integrated climb and descent segments. 

Computes a cruise flight path given the initial and final 
altitude and Mach number, initial weight, and desired 
range. 

Computes a reference flight path that follows an indicated 
airspeed (VIAPl) from sea level to 10000 ft, an indicated 
airspeed (VIAP2) from 10000 ft to Intersection with Mach 
number (RM3) , and RM3 up to cruise altitude; 

Sets up reference trajectory, either computed or read in. 

Computes coefficients for a continuous angle-of-attack 
perturbation command control law that maintains airspeed 
as a function of altitude. This is a closed-loop com- 
mand algorithm for climb ; 

Computes coefficients for a continuous angle-of-attack 
perturbation command control law that maintains altitude 
as a function of range-to-go. This is a closed-loop 
command algorithm for descent; 

Computes true airspeed, energy, energy rate, and fuel 
rate from indicated airspeed (or Mach number), altitude, 
and weight; 
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Under Category 2, the routines are: 


DATTRI 

DATTWN 

CDRAG 

CDRAG2 

CDRAG3 

CLIFTT 

CLIFT2 


CLIFT3 

ENGEPR 

ENGEP2 


ENGEP3 

ENGIDL 


FSUB 


TRIM 


Block data containing engine data for the tri-jet 
turbofan engine. 

Block data containing engine data for the twin-jet 
aircraft. 

Calls appropriate routine to compute the drag coefficient. 

Computes the drag coefficient for the tri-jet aircraft. 

Computes the drag coefficient for the twin-jet aircraft. 

Calls appropriate routine to compute the lift coefficient. 

Computes the lift coefficient as a function of Mach 
number, altitude, and angle-of-attack for the tri-jet 
aircraft. 

Computes the lift coefficient for the twin-jet aircraft. 

Calls appropriate routine to compute engine thrust and 
fuel flow rate. 

Computes the engine thrust and fuel flow rate as functions 
of altitude, Mach number, temperature variations, and EPR 
setting for the tri-jet aircraft. 

Computes the engine thrust and fuel flow rate for the twin- 
jet aircraft. 

Computes the engine thrust and fuel flow rate as functions 
of altitude and Mach number when EPR has been set at idle. 

Computes the derivative values of the first order differen- 
tial equations representing the longitudinal dynamics and 
fuel bum of the twin-jet aircraft. 

Computes the thrust and angle-of-attack for maintaining 
constant speed levels for a given altitude and cruise 
weight. 
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Under Category 3, the subroutines are: 

ATLOW Generates atmospheric density, pressure, temperature, and 
speed-of-sound as functions of altitude. 

WIND Computes the wind vector and its effect along the ground 

track of the aircraft. 

WlNDl Computes the longitudinal wind gradient as a function of 
altitude. 

WINDIN Reads in the data and sets up the wind profile as a 
function of altitude. 

WINDSH Computes wind gradient components as functions of altitude. 

Under Category 4, the subroutines are: 

DBLSRC Performs a linear double table look-up. 

FIAS Converts indicated airspeed in feet/second to Mach number. 

FIASM Converts Mach number to indicated airspeed in knots . 

GO The fourth-order Runge-Kutta-Gill numerical integration 

subroutines . 

OSUB Called by GO for special printout and to stop the integra- 

tion process when a variable reaches a certain magnitude. 

PAGE Starts a new page of printout. 

POLYEl Evaluates a polynomial for some fixed value of the 
independent variable. 

P0LY2 Evaluates a polynomial function of two independent variables . 

SERCHI Searches for a point in a monotonically increasing array. 

SGLSRC Performs a linear table look-up. 

TRACIT Traces subroutine calling sequence in case of program error. 

The interrelationship between these subroutines is also presented in 

Appendix B. 
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APPENDIX A 


AIRCRAFT EQUATIONS OF MOTION AND AUTOPILOT MODELS 


The objective of the TRAGEN program is to simulate an aircraft being 
steered to fly along an input or computed reference trajectory. This 
trajectory may be any combination of climb, cruise, and descent profiles. 
This simulation must be accurate enough such that the performance of the 
aircraft (in terms of fuel burned and time required to reach the destina- 
tion point) is adequately determined, as measured from the output. 

The purpose of this appendix is to provide the analytical expressions 
upon which the simulation was developed; this is done in three parts. The 
first section below defines the overall system and presents the differential 
equations of motion and fuel bum. The second section describes different 
methods for generating typical guidance commands and autopilot equations 
used for climb and descent. The third section derives the Breguet equa- 
tion used for cruise segment calculations. 


Equations of Motion and Fuel Burn 

To examine the vertical profile of the aircraft (i.e., altitude and 

airspeed vs range) , the longitudinal equations of motion are of primary 

importance. The control variables in this longitudinal plane are the 

angle-of-attack a and the magnitude of the thrust vector T. These quantities 

are shown with respect to aircraft airspeed V , lift L, drag D, weight 

s. 

W, and flight path angle y in Fig. A.l. 

The kinematic equations of motion of the aircraft in the longitudinal 
plane are 

X = Vg (A.l) 

h = V sin Y , 

a ' ’ 

where 

X = distance, or range, measured on the ground, 

h = altitude. 
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Figure A. 2. Elements of the Longitudinal Aircraft System 
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Vg = ground speed (magnitude of ~ ^ ^ ’ 

V = wind speed, 
w 


The inertial speed along the airspeed vector V is 

d 

V_ = V +V cos Y cos 5, (A. 2) 

law ’ ' 

where 6 is the angle between the horizontal projection of and V^. 

From Fig. A.l, the time rate of change of this vector for constant y is 

= V + V cos Y cos 6 = — (T cos a - D - W sin y) • (A. 3) 

law' m ' 


The time rate of change of the wind speed is 


V 

w 


-^h 

3h 

3V 


w 


3h a 


V sin y . 


(A.4) 


By substituting Eq. (A.4) into Eq. (A. 3) and solving for V , one obtains 

Qi 


• 1 

V = — (T cos a - D - W sin y) - rr^ V sin y cos y cos 6 . 

am ' 3h a ' ' 


(A.5) 


This Ignores the time rate of change of the wind direction. From Fig. A.l, 
one can write 


h = ^ (L cos y - W + T sin (y + a) - D sin y) . (A. 6) 

Equations (A.5) and (A. 6) represent the kinetic equations of motion of 
the aircraft. 


The remaining term that must be accounted for is the time-varying 
weight of the aircraft. Specifying the thrust also specifies the fuel 
burn rate f. Thus, the weight changes at the rate 

W = -f . (A. 7) 

Equations (A.l), (A.5), (A. 6), and (A. 7) are the five basic equations 
used to represent the longitudinal dynamics of the aircraft. 
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Further refinement could be added to these equations to include the 
effects of the following: 

1) . throttle dynamics (Including transient fuel flow rates) ; 

2) . relationship between throttle position, EPR setting, and thrust; 

3) . short period dynamics relating time rate of change of angle- 

of-attack, pitch rate, and pitch angle to elevator deflection; 

4) . required turning (lateral) motion for flying over fixed waypoints; 

and 

5) , variable wind direction and gust effects. 

However, these effects are considered to be of second order, and not re- 
quired for the intent of this simulation. For a more exact autopilot 
simulation, they would be required. 

The flight path angle is defined as 

Y = sin ^ (h/V ) (A. 8) 

3 

By differentating this expression and using Eqs (A. 5) and (A. 6) , one obtains 

1 2 

Y = — = — (T sin a - W cos y + L + V sin y cos 6 ). (A. 9) 

'mV ' 9h a 

a 

Equation (A. 9) can be used in place of Eq. (A. 6). 

Steering Procedures 

The climb and descent reference trajectories which are given to be 
followed consist of a sequence of points containing values of time, range, 
altitude, airspeed, flight path angle, specific energy, weight, and other 
variables. Any of these quantities which is measureable and monotonically 
changing can serve as the independent variable. To minimize airborne com- 
puter memory requirements, it is important to make the stored data represent- 
ing the reference trajectory as compact as possible. 

In this study, a set of steering equations are used to take points 
from the reference trajectory, convert these points to reference trajectory 
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commands, and then use these commands to set values of the control vari- 
ables. This steering process represents a rudimentary form of an autopilot. 

The steering process consists of commanding the thrust T and angle- 
of-attack a values so that the aircraft follows the reference as closely 
as possible. The system that Includes this process is depicted by the 
block diagram in Fig. A. 2. Note that flying along a reference trajectory 
consists of steering to connect a series of reference points. When a re- 
ference point is reached, new steering commands must be Issued so that the 
aircraft will then be guided to the next reference point. 

To fly along the reference path, an independent variable is first 
chosen. For this study, two different independent variables were chosen - 
range and altitude. Then, the remaining variables - primarily airspeed, 
flight path angle, and altitude (for range as the independent variable) - 
are stored as tabular functions of the chosen independent variable. 

Also, it is possible to fly along a nominal path using two approaches: 

1) . An open-loop approach where the thrust vector is directed 

over the next period in such a way that by the end of that 
period, the next reference point is reached. 

2) . A closed loop approach where the aircraft is continually 

steered to a continuously commanded trajectory which connects 
the reference points . 

Both of these approaches were examined for simulation of flying the climb 
profile. The closed loop approach gave superior performance, so only this 
approach has been retained. 

The problem with open loop steering was that it assumed that constant 
or linearly varying controls would cause the end points of a reference pro- 
file to be connected. This assumption did not account for perturbations 
due to wind, etc. along the way to be taken into account. Although the 
open loop methods produced paths which had roughly correct values of air- 
speed and altitude at given range values, these paths had large excursions 
from the reference flight path angle for the climb profiles. 
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Another problem with the open-loop approaches was that both a and T 
were varied to achieve fixed values of V and h for given range points. 

For optimum climb, thrust is usually set at the maximum value. Thus, usually 
only a remains as a valid control variable. 

Another consideration for implementing the climb profile is that there 
is no reason why a particular cruise condition (altitude, airspeed) has to 
be achieved when a certain range x is reached. Thus, a more logical Independent 
variable is altitude, with range allowed to be a free variable. 

For these reasons, a closed-loop steering approach was devised where 

reference values of flight path angle (with respect to the air mass) and 

airspeed are obtained as functions of altitude. (This assumes that altitude 

is monotonlcally Increasing during climb.) A perturbation control law was 

set up so that variations in a from a reference value were proportional to 

variations in Y and V from their respective command values. 

a 

For dynamic trim, when no wind shear is present, Eqs . (A. 5) and (A. 9) are 

mV = T cos a - D - W sin Y , (A. 10) 

mV y = T sin a + L - W cos Y = 0. 
a' ' 

That is, these non-linear equations must be continually solved for T and 
a to provide a dynamic condition where the specified acceleration is 
achieved for a steady flight path angle Y- 


Because y and V tend to change linearly with time, they can be considered 
Q. 

as ramp functions. Thus, the closed-loop controller should be considered to be 
at least a Type 1 system. From Eqs. (A. 10), the system perturbation equations 
are 

m 6V^ = - T sin a 6a - <5a - 1^ 6V^ - W cos ySy , (A. 11) 

m V 6 y = T cos a 6a + |^ <Sa + 6V + W sin y<5y • 

a ' 3a 3V a 
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The resulting transfer functions between y» V , and a are of the form 


6a 


(s/u)^ + 2^(s/w) + 1 



6a 


G (t s + 1) 

C C 2 

(s/oj)^ + 2c(s/oj) + 1 


(A. 12) 


where the time constants, etc. are functions of the parameters in Eq. (A. 11). 

The control problem can now be interpreted as shown in Fig. A. 3. To 
obtain the Type 1 system, the control law has to be of the form 

K K 

6a = (K +-h (V - V ) + (K +-^) (Y^ - Y), (A.13) 


where V and y are the commanded values of V and Y* This is the classical 
proportlonal-plus-lntegral controller. Gains are chosen to produce the 
desired response for removal of profile errors. 


For the climb profile, the subroutine STEERl mechanizes the above 

approach. To generate the continuous commands V and y , the computations 

ac c 

made at each reference point are 


3h 


3V 


3h 


\ + 1 - \ 
\+l " \ 


a , , - a 
n + 1 n_ 

h , - h 

n + 1 n 


Then , 



Y + (h — h ) 
'n n 



+ (h - h^) 


3h ’ 

3V 
a 

3h • 


(A. 14) 


(A. 15) 
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Figure A. 3 Control Loops for Perturbation Control of 
Airspeed and Flight Path Angle. 
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When the flight path angle Is very small (during the initial period of 
flight and when the aircraft levels off at 10000 ft to gain speed before 
resuming climb), Eqs. (A. 15) do not work well. For these cases, it is more 
appropriate to set 

0., (A. 16) 

V 

^ N+1 , 
and use the control law 

K, 

6 a = (K_ + -^) (Y - Y ) . (A.17) 

s c 

The above method, implemented as STEERl, worked quite well in causing the 
simulated profile to closely follow the reference path. Only one set of 
gain values was sufficient for the entire trajectory. 

For descending flight, the thrust again is usually constrained (idle) 
for optimum performance. Also, for this case, the main concern is to 
reach a fixed altitude when range-to-go to the destination point is a 
certain value. Thus, above 10000 ft, the airspeed can be allowed to be 
a free variable. For this case, only Inertial flight path Yj^ is required 
to be controlled. 



To generate a continuous command 
reference point is 


the computation made at each 


''ic ■ [<Vl - '■n> ' '’'n+l ■ ' 

Then, the control law is similar to Eq. (A.17), i.e. 

^4 

6 a = (K 3 + — ) (Yi^ - Yj) . 


(A. 18) 


(A. 19) 


where inertial values of flight path angle are used rather than those with 
respect to the air mass. Equations (A.18) and (A. 19) form the basis for 
the subroutine STEER2 which is used for closed-loop control of descending 
flight. Again, one set of gains is sufficient for the entire descent pro- 
file. 
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Cruise Computations Using the Breguet Equations 


A single cruise leg takes place in one vertical plane. Over this 
leg, it is assumed that the flight path angle is very small and that 
speed and altitude changes are negligible. Also, for now, it is assumed 
that there is no wind. With these assumptions specified, the following 
equations are valid. 

T = D , 

L = W , (A.20) 

X = V . 

a 

In addition, the time rate of change in weight can be expressed by the 
equation 

W = - T S , (A. 21) 

FC 

where is the engine specific fuel consumption, 
r L 


These equations can be used to formulate the standard range equation 
as follows; 

X = (dx/dW)W = V (A. 22) 


Therefore, 

dx/dW 



-V 


T(s"^J 


-V^(L/W) -V^(L/D) ^ 

(T/D)D(Sp^) " W 


(A. 23) 


The Breguet factor or range factor, R , is defined as: 

r 




V^(L/D) 


Then, 


dx/dW = -Rp ^ . 


(A. 24) 


(A. 25) 
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From Eq. (A. 25), one can write 


X = 


i 


^final ^ 


W 


dW 


initial 


or 


— „ / ^initial \ 




(A. 26) 


A(.27) 


where Rp is the average value of over the range traveled. Using the 
average value for the range factor Rp is an approximation but a very good 
one for cruise performance. Equation (A. 27) is referred to as the range 
equation. 


The range equation is often used to determine an optimum altitude and 
Mach number to maximize the range. However, for the purpose of the TRAGEN 

program, cruise speed, altitude and the required range of the cruise seg- 
ment are specified, and it remains to find the fuel burn over the cruise 
segment. Thus the range equation is rewritten as follows: 


First, the fuel burned is 

Fuel initial final 

Then, the range traveled is 


(A. 28) 


X 


Rp In 


1 - W. t/W. , 

fuel initial 


(A. 29) 


Thus, the fuel burned to achieve a given range x is 

"fuel = (A.30) 

* Note that this is a relatively trivial optimization result for a 
commercial transport aircraft since the cost of time is not con- 
sidered and the climb and descent legs are ignored in the problem. 
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In TRAGEN, the average Breguet factor is computed by evaluating 
Eq. (A. 24) at the initial and final altitude and airspeed conditions 
specified to be achieved over the given range. Equation (A. 30) is used 
to iterate on the amount of fuel burned over this segment. This is used 
in turn to compute the final weight to determine the trimmed value of 
Eq. (A. 20) and to obtain the final value of R^. 
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APPENDIX B 


TRAGEN SUBROUTINE DESCRIPTION 


This appendix contains an explanation of the data storage for program 
OPTIM. Following that is an explanation of the main program, the sub- 
routines, and the functions in alphabetical order. 

Data Storage 

The major part of the data communications between subroutines in 
OPTIM is through labelled common statements. There are ten such commons. 
Their names and a short description of each are: 

CODE Cruise, climb, descent variables. 

CONST Constants. 

CRUISE Cruise table and associated variables. 

DESCRP Assorted variables. 

ERROR Error traceback information. 

GRAPH Data to be written to Unit 11 and associated variables. 

Includes the final climb and descent trajectories. 

INPUT All input parameters. 

TOA Time-of- arrival and step climb variables. 

TRIJET Engine data, tri-jet. 

TWINJT Engine data, twin-jet. 

WINDP Wind input data and associated parameters. 

As a convenience, the CDC UPDATE capability is used to insert COMMON 
statements into source decks. This facilitates changing items in COMMON 
with no loss of program portability, because UPDATE produces a compile file 
which is directly llstable, editable, and compilable by any standard 
FORTRAN. 

Program Explanation 

Following is a description of the TRAGEN routines. 
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MAIN PROGRAM; TRAGEN 


The sequence of steps in TRAGEN follows the flow chart presented in 
Section IV. The purpose of TRAGEN is to simulate an aircraft being com- 
manded to follow a reference profile or to cruise. The profile may either 
be a climb or a descent as governed by the mission segment flag NSC. It 
may be either input or computed as governed by the flag ITRAJ . 


The program first reads in control flags and other data. This 
is followed by reading in of variables to initialize the aircraft state. 

If wind is to be used, it is read in before the initial segment is computed. 
Likewise, if a reference trajectory is to be read in from Unit 11 it is 
brought in for the Initial segment. 

If the segment is an initial climb, or a descent, the control law 
gains are read in. Then subroutine SETREF is called to set up the 
reference. CMACTL is then called to integrate the trajectory for com- 
parison with the reference. 


If the segment is a climb following some other segment , a step climb 
to the desired final altitude is inserted first. TRAGEN then continues 
with the climb and descent logic. 

If the segment is a cruise, subroutine ACRUSE is called to control 
the computation. 


After each segment, TRAGEN returns to the beginning to read a new 
header and new options (as a minimum) . It must be ended with a header 
card followed by an NSC = 5 card. 


Subroutines called: 


Commons used: 


ACRUSE 

ATLOW 

CMACTL 

SETREF 

TRACIT 

TRIM 

WINDIN 

WINDSH 


CCDE 

ERROR 

INPUT 

TQHCOM 

TRAGEY 

TVAR 
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ACRUSE 


Subroutine ACRUSE prints the initial output for the cruise segment, 
sets up the initial conditions depending on whether or not this is the 
initial segment, and calls subroutine CREWZ to perform the cruise com- 
putations. If IWIND is not zero, ACRUSE reads in the cruise wind. 


Subroutines called: 


Commons used: 


ATLOW 

CRUSE 

FIAS 


CCDE 

ERROR 

INPUT 

TRACEY 

TVAR 

WINDP 


ATLOW 

This subroutine generates the atmospheric density (in lb sec^/ ft^) , 
atmospheric pressure (in Ib/ft ) , atmospheric temperature (in degrees 
Kalvin) and speed of sound (in ft/sec) at a given altitude below 20,000 
meters (65,617 feet). It also makes the appropriate modifications in 
atmospheric density and speed-of-sound to account for variations in 
standard day temperature (represented by the input DTEMP). The 1962 
standard atmosphere is used. This version of the program does not 
calculate a new atmosphere when called at successive times at the same 
altitude. 


CONDATA 

This subroutine contains all the data for program constants. 
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BLOCK DATA - DATTRI 


This block data contains the engine data used with the tri-jet aircraft 
model. Three tables are used to describe idle thrust, idle fuel flow, and 
maximum continuous engine pressure ratio (EPR) . 

Subroutines called: Common used: 

None TRIJET 


BLOCK DATA - DATTWN 

This block data contains numerical characteristics of the turbofan 
engine used with the twin-jet aircraft model. Seven tables are used to 
describe idle thrust and fuel flow for bleed valves open and closed, 
altitude of surge bleed valve closure, maximum EPR for climb and cruise, 
and Mach number corrections. 

Subroutines called: Common used: 

None TWINJT 


CDRAG 


This subroutine calls the appropriate routine to compute the aircraft 
drag coefficient CD based on the particular aircraft model selected by 
the input variable lAC. Currently, two models are available, but logic 
is present to use up to four different aircraft. 

Subroutines called: Commons used: 

CDRAGl* None 

CDRAG2 

CDRAG3 

CDRAG4* 

* not included with program 
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CDRAG2 


This subroutine computes the drag coefficient CD for some given 
Mach number EM and lift coefficient CL for a medium range tri-jet trans- 
port aircraft model. The value is computed from the coefficients of a 
polynomial stored in the array COEFF. 

Subroutines called: Commons used: 

P0LY2 None 


CDRAG3 

This subroutine computes the drag coefficient CD for some given 
Mach number EM and lift coefficient CL for a medium range twin jet trans- 
port aircraft model. CD is computed by polynomial evaluation, including 
interpolation of the polynomial and its first derivative in certain 
regions, as necessary. 

Subroutines called: Commons used: 

P0LY2 None 


CLIFTT 


CLIFTT calls the appropriate routine to compute the aircraft lift 
coefficient for the particular aircraft model selected by the input 
variable lAC. Currently, two models are available, but logic is present 
to use up to four different aircraft. 

Subroutines called: Commons used: 

CLIFTl* None 

CLIFT2 

CLIFT3 

CLIFT4* 

* not included with program. 
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CLIFT2 


This subroutine computes the lift coefficient CL for a medium range 
tri-jet transport aircraft as a function of Mach number EM, altitude H, 
and angle-of-attack ALPHAP. The lift coefficient consists of three terms: 

Cl = (basic) + + Cj^ a 

The first term C (basic) is a polynomial function of angle-of-attack a. 

L« 

The value of this term is checked against the buffet boundary expressed 
as a polynomial of Mach number. The second term C^^^ is a polynomial of 
Mach number with altitude as the parameter. The third term C^^ is also 
a polynomial of Mach number. The coefficients of the polynomial are fit 
for different altitudes. 

Subroutines called: Commons used: 

POLYEl None 


CLIFT3 

This subroutine computes the lift coefficient CL for a medium range 
twin jet transport aircraft as a function of Mach number EM, altitude H, 
and angle-of-attack ALPHA. The lift coefficient consists of three terms: 

Cl = C^ (basic) + + C^^ cx 

The first term C (basic) is a function of angle-of-attack. The second 

term C, is a function of altitude and Mach number. The third term C, 

Lo La 

is also a function of altitude and Mach number. These terms are determined 
by table lookup. 

Subroutines called: 

DBLSRC 
SERCHI 


Commons used: 
None 
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CMACTL 


This subroutine computes the actual trajectory for comparison with the 
reference trajectory. First, CMACTL initializes all variables that are 
printed out and that are modified by the integration process. At this 
point, the update process is ready to begin. 

The first step of the update (integration) process is writing the 
simulated, time-varying state variables as determined from integrating 
the equations of motion. These are followed by a written line of variables 
obtained from the reference trajectory at about the same point along the 
profile. The reference trajectory data points are separated by steps in 
specific energy or altitude of 500 ft. 

The next step is to determine whether the end of the integration 
process has been reached. The subroutine will exit when any of the following 
take place : 

time T ^ TSTOP, 

counter ICT > NOPT, 

where NOPT is the number of points in the reference trajectory. 

Next, the process of generating the steering commands to follow the 
reference trajectory is simulated. One option (STEERl) is present for 
climb commands and one (STEER2) is present for descent command generation. 
These commands consist either of airspeed and/or flight path angle which 
are used to command continuous changes to angle-of- attack. They also 
compute how long (DT) these commands hold until the next set of commands 
should be issued. 

Following the issue of the steering commands, the aircraft equations 
of motion are integrated by calling the integration subroutine GO. Then 
the counter ICT is updated, certain output variables are computed, and 
the program loops back to begin the update cycle again. 


47 


Subroutines called: 


Commons used: 


ATLOW 

GO 

STEERl 

STEER2 


CCDE 

ERROR 

INPUT 

TQHCOM 

TRAGEN 

TVAR 


CREWZ 

Subroutine CREWZ calculates the cruise performance of the aircraft 
when initial and final altitude and speed, initial weight, and desired 
range are given. An iteration employing Breguet factors is used. 
Appendix A describes the theory upon which this subroutine is based. 


Subroutines called: 


Commons used : 


ATLOW 

FIASM 

PAGE 

TRIM 

WIND 


CCDE 

ERROR 

INPUT 

TRAGEY 


DBLSRC 

This function performs a double table lookup. Given a function f(x,y), 
this function interpolates the appropriate arrays to obtain approximate 
values of f(A,B). The four points which surround (A,B) are first found, 
and the function is evaluated at each. Then these values are interpolated, 
first on x and then on y, to obtain the approximate solution. 

Subroutines called: Commons used: 

SERCHI None 
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ENGEPR 


This subroutine calls the appropriate routine to compute the aircraft 
maximum thrust and EPR, the thrust associated with the input EPR, and the 
fuel flow rate. The engine model is associated with the particular air- 
craft model selected by the input variable lAC. Currently, two models 
are available, but logic is presented to use up to four different aircraft. 

Subroutines called: Commons used: 

ENGEPl* None 

ENGEP2 

ENGEP3 

ENGEP4* 

* not included in program 
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ENGEP2 


This subroutine generates the thrust THRUST and the fuel flow rate 
FDOT for some given altitude H, Mach number EMAKNO and EPR setting. First 
EPRMX, the maximum continuous EPR, is determined by table look-up 
for some given temperature Ta and altitude H, where 

Ta = T(1 + (EMAKNO)^)^ - 273.15. 

Here, T is the temperature corresponding to altitude H, after temperature 
variation correction, and 

Y = 1.4, the ratio of specific heats. 

The EPR setting is limited to EPR ^ EPRMX for cruise and EPR ^ EPRMX - . 1 
for climb or descent. 

Second, (FN/6^) is computed from a polynomial. Then, the thrust is 
computed as, 

THRST =3(6 ) (FN/6 ). 

am e 

This is the thrust for the medium range tri-jet transport aircraft model. 

Here, 6 is the pressure ratio 
am 



o 


Here, P is the atmospheric pressure corresponding to altitude H, and P^ is 
the atmospheric pressure at sea level. A factor of 3 is used since there 
are three engines. 

Finally, the fuel flow rate FDOT is computed as: 


FDOT = 3 *WFC * 6a *Kc 

where 

Kc = ,00 223181 Ta + .9675897, 

5a = 6 (1 + (EMAKNO) 

am z 
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Also, WFC is the fuel-flow rate computed as a polynomial of EPR, where the 

coefficients of the polynomial depend on both altitude and Mach number. 

Subroutines called: Commons used: 

ATLOW CCDE 

DBLSRC ERROR 

POLYEl INPUT 

TRIJET 
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ENGEP3 


This subroutine generates the thrust THRUST and the fuel flow rate 
FDOT for some given altitude H, Mach number EM and EPR setting. First 
EPRMX, the maximum continuous EPR, is determined by table look-up, 

(Tables 6 and 7 in Block Data) for some given temperature Ta and altitude 
H, where 

Ta = T(1 (EM)^) - 273.15. 

Here, T is the temperature corresponding to altitude H, after temperature 
variation correction, and 

Y = 1.4, the ratio of specific heats. 

The EPR setting is limited to EPR ^ EPRMC. 

Second, (FN/6^) is computed from a polynomial. Then, the thrust is 
computed as, 

THRST =2(5 ) (FN/6 ) , 

am e 

where 6 is the pressure ratio 
am 



o 


Here, P is the atmospheric pressure corresponding to altitude H, and P^ is 
the atmospheric pressure at sea level. A factor of 2 is used since there 
are two engines. 

Finally, the fuel flow rate FDOT is computed. A polynomial is used 
to calculate WFC for a given EPR and altitude. At values of EPR < 1.6, 
there is also a correction for Mach number (Table 10 in Block Data) . 

Then: 

FDOT = 2 *WFC * 6a *Kc , 

where 

K = .0022 T + 0.97 , 

c a ’ 

Sa . 6 (1 + . 

am 2 
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Subroutines called: 


Commons used: 


DBLSRC 
POLY 2 
SGLSRC 


CCDE 

ERROR 

INPUT 

TWINJT 


ENGIDL 


ENGIDL is called during descent to compute thrust and fuel flow 
rate for idle EPR. It does this through table look-up for the appro- 
priate aircraft. 


Subroutines called: 


Commons used: 


DBLSRC ERROR 

SGLSRC TRIJET 

TWINJT 


FI AS 

FIAS returns Mach number as a function of indicated airspeed (in 
feet per second) and atmospheric pressure. 

Subroutines called: Commons used: 

None None 


FIASM 


FIASM returns indicated airspeed in knots as a function of Mach 
number and atmospheric pressure. 


Subroutines called: 


Commons used: 


None 


CONST 
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FSUB 


This subroutine computes the derivatives for the Integration process. 
It Is called four times per Integration step by the Integration subroutine 
GO. The following steps are taken to compute the derivatives: 

1. The flight path angle with respect to the air mass Is computed as 

Y = sin ^ h/V 

2. The Mach number and atmospheric density are computed by calling 
ATLOW. 

3. The longitudinal wind magnitude and Its gradient with respect to 
altitude are computed by calling WlNDl. 

4. The control variables are computed based on the steering option 
used. They are for ICNFL: 


STEERl option. Here, airspeed and fllght-path-angle commands/ 
are computed as linear functions of altitude: 


^c ” ^cl ^ 3h ’ 

+ "Ih- 


Then, the actual values are subtracted from these commands 
to generate 6y and 6V errors. A proportional plus Integral 
control law of the form 


6a = 5V + K 2 /6V + 6y + / <5y. 


Is used to compute continuous perturbations to the steady 
value of angle-of-attack. Angle-of-attack Is limited be- 
tween (-4°, + 22°). Thrust Is set at a maximum value. 
During the level parts of the climb (y ^ 0.1°), only the 
flight path errors (6y, /■5Y) are used In the control law. 
Each segment of this option Is usually cut off on time. 
However, during level acceleration, the segment Is cut 
when airspeed reaches the next reference value. 
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2: STEER2 option. Here, a constant inertial flight path 

angle (GMC) is commanded. The actual value is estimated 
and subtracted, and a proportional plus Integral control 
law is used to compute perturbations to the angle-of- 
attack command: 

6a = 6 y + / 6y. 

Again, total angle-of-attack is limited to (-4°, + 22°). 
The thrust is set at idle value by setting EPR to 1.1. 

5. Lift, drag, and mass of the aircraft are computed. 

6. Thrust and fuel flow rate are computed. 

7 . The five basic derivatives representing the longitudinal 
dynamics of the aircraft are computed: 


h = (L cos Y “ W + T sin (y + a) - D sin y)/m, 

3V^ 

V = (T cos a - D - W sin y) /m V sin y cos y cos 6, 

X = V 

g 

• • 

h = h , 

W = -f/3600. 


Subroutines called: 


Commons used: 


ATLOW 

CDRAG 

CLIFTT 

ENGEPR 

ENGIDL 

WIND 

WINDl 


CODE 

ERROR 

INPUT 

TRAGEY 

TQHCOM 

WINDP 
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GO 


This subroutine is a Runge-Kutta-Gill fourth order numerical integra- 
tion package which integrates a set of eight first order ordinary differen- 
tial equations. The step size of the independent variable is H. X and 
XF are the initial and final values of the Independent variable (which is 
time , for this application) . 

Subroutines called: Commons used: 

FSUB ERROR 

OSUB TQHCOM 


OSUB 


This subroutine has two purposes: (1) to write out intermediate values 
of system variables at the end of each integration step, and (2) to stop 
integration along a certain segment. For ICNFL: 

;1 STEERl option. For climbing flight (y > 0.1°; ICFL =1), the 
segment is cut off (XF = XDQ(l)) when altitude H reaches the 
next reference value HP. For near level flight (y ^ 0.1°; 

ICFL = 2) , the segment is cut off when airspeed VA reaches 
the next reference value VCl. 

:£ STEER2 option. For descending flight, the segment is cut off 
when the range value X becomes greater than the next reference 
value RP. 


Subroutines called: 
NONE 


Commons used: 

INPUT 

TRACEY 

TQHCOM 
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PAGE 


This subroutine advances the printout to the top of the next page. 


POLYEl 


This function evaluates the polynomial 

Y = b(l) + b(2) X + b(3)x^ = ...b(M)x““^ 


P0LY2 


P0LY2 evaluates the polynomial 


Z + . . . 


n-1 


+ C2j_Xj_ + C22Xj_X2 + . . . C2^Xj_X2 


n-1 


+ ... 


, in-1 , m-1 , . m-1 n-1 

+ c -X, + c oX, x„ + ... + c X- Xo 

ml 1 m2 1 2 mn 1 2 
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REFCOM 


This subroutine computes a reference flight profile that is similar 
to one that would be specified in a pilot's handbook for a particular 
aircraft. (For example, aircraft are usually limited to be under 250 kt 
IAS below 10000 ft. The tri-jet might have a climb schedule of 320 kt 
IAS/0. 73M and a descent schedule of 0.73M/320 kt IAS.) Thus, this subrou- 
tine computes a climb profile that follows the following sequence: 

1. Accelerate from VO to VIAPl (indicated airspeed in kt) at 
altitude HO. 

2. Climb to 10000 ft at VIAPl in 500 ft steps. 

3. Accelerate from VIAPl to VIAP2 (indicated airspeed in kt) 
at 10000 ft. 

4. Climb to intersection with Mach number RM3 at indicated air- 
speed VIAP2, in 500 ft steps. 

5. Climb to altitude HF at Mach number RM3, in 500 ft steps. 

At each step, the variables time (Tl) , range (Rl), altitude (HI), true 
airspeed (VTl), flight path angle (GAMl) , specific energy (El), fuel 
burned (FI), power (EPR) setting (EPl) , and windspeed are computed and 
stored in the array A. 

If the path is a climb profile, maximum EPR is used. For descent, 

EPR is set to 1.1 and xdle thrust and fuel flow rates are used. 

REFCOM calls the subroutine VTCM to convert indicated airspeed or 
Mach to True airspeed and to obtain thrust, energy, energy rate and other 
variables at a particular altitude. Then the computation sequence is 

AE = E - E 

P 

At = E/E 

Y = sin ^( (Ah/At) /V.J. ) 

P 

AR = (V + V ) At/2 
P 

W = W - W At 
P P 
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where the subscript p indicates the value of a variable at the previous 
altitude. 


Subroutines called: 


Commons used: 


VTCM CCDE 

WIND ERROR 

INPUT 

TRAGEY 

TVAR 

WINDP 


SERCHI 


The array TX(*) is monotonically increasing. This subroutine searches 
the index Z such that 


TX^ < X < TX^^^ . 


and returns both S. and pf where 
X - TX„ 


pf = 


^^£+1 - 
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SETREF 


This subroutine sets up the reference trajectory, either by reading 
it in or by calling REFCOM to compute it. If input parameter ITRAJ = 1, 
the data input procedure is set to accept output from the companion 
program OPTIM. For a climb profile, this comes in the form: 

CGRAF(I,J) for 1=1, JCLIMB and J = 1,10. 

JCLIMB is the number of data points. There are up to 10 variables for 
each point. For the descent profile, the input data come in the form: 

DGRAF(I,J) for 1=1, JDESCN and J = 1,10. 

JDESCN is the number of data points, and again, there are up to 10 
variables for each point. Because OPTIM generates the descent profile 
backwards in time, the DGRAF array variables are reordered with time and 
range given negative values, and fuel burned is manipulated to be sub- 
tracted from the initial weight rather than added to the final weight. 

If ITRAJ = 2, SETREF sets up the input and calls the subroutine 
REFCOM to compute the reference profile. Thxs may either be a climb or 
descent profile as governed by the flag NSC. 

Subroutines called: Commons used: 

PAGE ERROR 

REFCOM INPUT 

TRAGEY 

IVAR 


SGLSRC 


This function evaluates a single function F at the point A. This 
is done by linear interpolation to obtain A' s location in the array X 
and using the tabulated values of F(X) . 

Subroutines called: Commons used: 

SERCHI None 
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STEERl 


This subroutine provides air-referenced flight-path-angle and airspeed 
conmands that are used in FSUB for closed-loop steering during climbing 
flight. This routine is based on the assumption that thrust is set to 
maximum value and that angle-of-attack perturbation commands can be re- 
lated to the difference between actual and commanded values of airspeed 
and flight path angle. The flight path and airspeed commands are generated 
as functions of altitude from values (Y.»V,,h.,Y,V,h) taken from 
the reference trajectory: 


3h h , - h 
+ n 


9h 


V . - V 
+ n 

h^ - h 
+ n 


Y 


cl 




cl 


V -h IJ. 

n n 3h 




At = 2(h^ - h)/(V^ sin Y+ + sin y) 


When the next reference value of flight path angle (Y^) is less than 0.1°, 
the above equations are replaced with 



Subroutines called: 
NONE 


Commons used: 
TVAR 
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STEER2 


This subroutine provides inertial flight-path-angle commands that are 
used in FSUB for closed-loop steering during descending flight. This routine 
is based on the assumption that thrust is set to idle value and that angle- 
of-attack perturbation commands can be related to the difference between 
actual and commanded values of flight path angle. The flight path angle 
commands are generated by keeping altitude as a fixed function of range-to- 
go to the landing point. Values of altitude (h^) and range-to-go (r^) are 
taken from the next reference point. Then 


3h 


h 

3X 

r - 
+ 

X ’ 

Yc = 

tan 

1 .3h 
^3X 


This value of flight path angle command is stopped when the next reference 
value of range-to-go (r ) is reached. 

Subroutines called; Commons used: 

NONE TVAR 


TRACIT 

In case of error, this subroutine provides a "walk back" through 
the subroutine calling hierarchy. If the subroutine is set up to 
recognize the computation or logic to be in error, then TRACIT can 
be used to find the source of the error. 
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TRIM 


This subroutine is used to compute the trim conditions for 
medium renge transport aircraft. This subroutine computes angle~of~ 

attack a and thrust T to keep the aircraft in trim for constant speed 
level flight, for a given altitude and for a given Mach number. 

With Y the flight path angle, the equations of motion in the horizontal 
and vertical directions are as follows; 

— (dv/dt) = T cos a - D - W sin Y 
6 

— v(dY/dt) = T sin a + L - W cos y 
S 

For a trimmed condition: 

(dv/dt) = (dY/dt) = 0. 

The two equations are combined to eliminate thrust to give the equation: 
(WcosY - L)cosa - (Dsina + WsinY)slna = 0. 

This equation is solved by iterating with angle-of-attack, a. 

Once the aircraft is trimmed, the thrust is solved from the drag by 
T = (D + WslnY)/cosa . 

This required thrust is matched by iterating on values of power setting 
(EPR) and calling subroutine ENGEPR. Once the correct power setting is 
determined, the engine fuel flow is also known. 

Subroutines called; Commons used: 

CDRAG CODE 

CLIFTT DESCRP 

ENGEPR ERROR 

INPUT 


63 



VTCM 


This subroutine computes true airspeed, energy, energy rate, thrust, 
drag, lift, and fuel rate from indicated airspeed, altitude, and weight. 
The computation sequence is as follows: 


Pressure 
Temperature 
Density 
Speed of sound 

Mach number 


p = f(h) 

T = f(h) 

p = p/(3092.4 T) 
a = 65.76 /r 

M = (5.(((((V^g/2496.5)^ + l.)^'^-l,) 

(2116.22/p) +1.)^'^^ -1.))^^^ 


True airspeed - M a 

If M is given, then 


Indicated airspeed “ 2496. 5((( (p/2116. 22) 

((1. + .2M^)^'^-1.) +1.)^^^-!.))^^^ 


Specific energy 
Thrust 

Fuel flow rate 
Lift 


Drag 

Energy rate 


E = h + /2.g 

Th = f(h,M,EPR) 
w = f(Th,h,M) 

L = W 

= L/( S/2) 

S = 

D = pV^^ S Cp/2 
E = (Th - D)V^/W 


Subroutines called; 


Commons used: 


ATLOW FIAS 

CDRAG FIASM 

ENGEPR 
ENGIDL 


CCDE 

ERROR 

INPUT 
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WINDIN 


This subroutine rends in the wind profile (the magnitude and direction 
of wind as a function of altitude) , The wind magnitude input is in knots 
and the program converts it to ft/sec and stores it in the VW array. The 
wind direction is stored in PSIW in degrees. The input represents the 
direction the wind is coming from. The program adds 180“ to this value 
to obtain the vector direction. The altitudes corresponding to these 
wind magnitudes and directions are stored in array HWIND. The wind may 
be input as a single profile, valid over the entire flight, or as separate 
climb, cruise and descent profiles. In the case of a step climb, the cruise 
profile is used for lower cruise, step climb, and upper cruise segments. 

Subroutines called: Commons used: 

None INPUT 

WINDP 


WINDSH 


This subroutine uses the input profile of wind magnitude VW( ) and 
direction PSIW( ) to compute north and east components of wind shear as 
a function of altitude HWIND( ) . These shear components change every 
2000 ft of altitude. 

Subroutines called: Commons used: 

NONE INPUT 

WINDP 
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WIND 


This subroutine computes the wind velocity as a function of altitude. 
This is combined with the aircraft velocity with respect to the air mass 
to compute ground velocity. Inputs to this program are H, the altitude 
in feet; PSIG, the aircraft heading in degrees; VTAS, the aircraft air 
speed; GAMMR, the angle of attack; VW, and PSIW arrays. The outputs 
from this program are VWA, the wind speed, and VG, the aircraft ground 
speed. 

Subroutines called: Commons used: 

SERCHI DESCRP 

ERROR 

INPUT 

WINDP 


WINDl 


This subroutine computes the gradient of the wind speed in the north 
and east directions from the input array WINCP( ) and the altitude. These 
values are then used to compute the gradient of wind DWDH with respect to 
altitude along the longitudinal axis of the aircraft. The wind magnitude 
is also computed. 

Subroutines called: Commons used: 

INPUT 

SERCHI WINDP 


66 



1 Report No 2 Government Accession No. 

NASA CR-166062 

3 Recipient's Catalog No 

4 Title and Subtitle 

TRAGEN - COMPUTER PROGRAM TO SIMULATE AN AIRCRAFT 
STEERED TO FOLLOW A SPECIFIED VERTICAL PATH - 
USER'S GUIDE 

5 Report Date 
May 1983 

6 Performing Organization Code 

7 Author(s) 
Anon. 

8 Performing Organization Report No 

10 Work Unit No 

9 Performing Organization Name and Address 

Analytical Mechanics Associates, Inc. 
2483 Old Middlefield Way 
Mountain View, California 94043 

1 1 Contract or Grant No 

NASl-15497 

13 Type of Report and Period Covered 
Contractor report 

12 Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 

14 Sponsoring Agency Code 

15 Supplementary Notes 

Langley technical monitor: Kathy H. Samms 

Program User's Guide 

16 Abstract 


This document is a technical description and a user's guide for a computer 
program — called TRAGEN — which is used to simulate an aircraft steered to 
follow a climbing, cruise, or descending profile or any sequential combination 
of these flight phases. Specifically, the program simulates the longitudinal 
dynamics of a medium range twin-jet or tri-jet transport aircraft. For the 
climbing trajectory, the thrust is constrained to maximum value, and for descent, 
the thrust is set at idle. For cruise, the aircraft is held in the trim condition. 
For climb or descent, the aircraft is steered to follow either (a) a fixed profile 
which is input to the program or (b) a profile computed at the beginning of that 
segment of the run. For climb, the aircraft is steered to maintain the given 
airspeed as a function of altitude. For descent, the aircraft is steered to 
maintain the given altitude as a function of range-to-go. In both cases, the 
control variable is angle-of-attack. The given output trajectory is presented 
and compared with the input trajectory. Step climb is treated just as climb. 

For cruise, the Breguet equations are used to compute the fuel burned to achieve 
a given range and to connect given initial and final values of altitude and Mach 
number. 


17 Key Words (Suggested by Author (s)) 
Traiectory simulation 
Fuel conservation 
Flight management 

18 Distribution Statement 

Unclassified - 

- Unlimited 

19 Security Ctassif (of this report) 
Unclassified 

20 Security Classif (of this page) 
Unclassified 

21 No of Pages 

67 

22 Price 


11-305 


For sale by the National Technical Intormation Service, Springfield, Virginia 22161 




















End of Document 



